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A phenomenological theory of magnetic states in noncentrosymmetric tetragonal antiferromagnets 
is developed, which has to include homogeneous and inhomogeneous terms (Lifshitz-invariants) 
derived from Dzyaloshinskii-Moriya couplings. Magnetic properties of this class of antiferromagnets 
with low crystal symmetry are discussed in relation to its first known members, the recently detected 
compounds Ba2CuGe207 and K2V3O8. Crystallographic symmetry and magnetic ordering in these 
systems allow the simultaneous occurrence of chiral inhomogeneous magnetic structures and weak 
ferromagnetism. New types of incommensurate magnetic structures are possible, namely, chiral 
helices with rotation of staggered magnetization and oscillations of the total magnetization. Field- 
induced reorientation transitions into modulated states have been studied and corresponding phase 
diagrams are constructed. Structures of magnetic defects (domain-walls and vortices) are discussed. 
In particular, vortices, i.e. localized non-singular line defects, are stabilized by the inhomogeneous 
Dzyaloshinskii-Moriya interactions in uniaxial noncentrosymmetric antiferromagnets. 
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I. INTRODUCTION 

In many magnetic crystals the magnetic properties 
are strongly influenced by the antisymmetric exchange 
{Dzyaloshinskii-Moriya) coupling which is generally de- 
scribed by a vector product formed by the magnetic mo- 
ments Si of two magnetic ions: 



wd = Djj • (Sj x Sj) 



(1) 



and the so-called Dzyaloshinskii vector D^oB Based on 
phenomenological considerations, the interaction ([!]) was 
introduced by Dzyaloshinskii to explain the observation 
of a small net magnetization in a number of antifepromag- 
nets, a phenomenon called weak ferromagnetisma which 
is due to a slight deviation of the sublattice magnetiza- 
tions from antiparallel arrangement. Extending Ander- 
son's theory of superexchange Moriya later found a mi- 
croscopic mechanism due to spin-orbit interactions that 
is responsible for the interactions (^). They arise in 
certain groups of magnetic crystals with low symmetry 
where the effects the couplings ([[]) do not cancel.El Dur- 
ing the following decades intensive theoretical and experi- 
mental studies on the Dzyaloshinskii-Moriya coupling (|l|) 
resulted in a deep insight into its microscopic origins and 
its manifestation in macroscopic properties of magnetic 
materials. oH Now it is known that weak ferromagnetism 
must essentially be attributed to many types of complex 
magnetic structures. It influences appreciably the mag- 
netic properties of several important classes of magnetic 
materials such as orthoferrites, manganites, some bifih- 
temperature superconducting cuprates, and others.Qu 

Another fundamental macroscopic manifestation of an- 
tisymmetric couplings (|l|) takes place in noncentrosym- 
metric magnetic crystals. Dzyaloshinskii showed that, in 
this case, the interaction (Q) stabilizes long-periodic spa- 
tially modulated structures with fixed sense of rotation 



of the vectors Within a continuum approximation 
for magnetic properties, the interactions responsible for 
these modulations are expressed by inhomogeneous in- 
variants. We will call these contributions to the (free) 
magnetic energy, involving first derivatives of magnetiza- 
tion or staggered magnetization with respect to spatial 
coordinates, inhomogeneous Dzyaloshinskii-Moriya inter- 
actions. They are linear with respect to the first spatial 
derivatives of the magnetization M of typeO 



*, dM 3 

an 



M, 



dMj 
drj 



(2) 



where Mi, Mj are components of magnetization vec- 
tor(s) that arise in certain combinations in expressions 
(^) depending |-pn crystal symmetry, and 77 is a spa- 
tial coordinated Such antisymmetric mathematical forms 
first have been studied in the theory of phase transitions 
by E. M. LifshitzEl and are known as Lifshitz invariants. 
The first magnetic modulated chiral structures predicted 
inl3 were observed in the cubic noncentrosioametric crys- 
tals MnSi and FeGe (space group P2i3).Qo During the 
following years modulated magnetic structures of this 
kind were discovered and investigated in several.plasses of 
magnetic crystals lacking inversion center.OEJiiij These 
chiral helical structures are essentially different from nu- 
merous other spatially modulated magnetic states in sys- 
tems with competipgpexchange interactions (as, e.g., in 
rare earth metals) Ho The latter are characterized by 
rather short periods (usually including only few unit 
cells) and arbitrary rotation sense. On the contrary, chi- 
ral structures due to ([!]) or (j^) have long period and a 
fixed sense of rotation. For example, in MnSi the period- 
icity lengths of the helix in zero magnetic field was found 
to be about 170 A (39 unit cells), andpFeGe has an even 
larger period (700 A or 149 unit cells). □ The interactions 
of type (0) may also stabilize periodic structures mod- 
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ulated in two-dimensions (uortex lattices) and localized 
axisymmetric inhomogeneous states.t£l 

Up to now both physical effects induced by 
the antisymmetric exchange coupling ([l]) ■ weak- 
ferromagnetism and chiral modulations — never have 
been observed simultaneously in one magnetic system. 
Moreover, in noncentrosymmetric magnetic crystals with 
chiral modulations that were known and described so 
far, the existence of weak ferromagnetism is excluded 
because of their symmetry.tll In this paper, we show 
that both phenomena can coexist in the recently dis- 
covered nonceiitaisymmetric tetipgonal antiferromagnets 
Ba 2 CuGe 2 7 No and K 2 V 3 8 M Due to the crystallo- 
graphic and magnetic structures of these crystals, the 
Dzyaloshinskii-Moriya coupling (|l|) favours noncollinear- 
ity along one direction and spatial modulations along the 
others. Here we determine possible magnetic phases and 
study their evolution in applied magnetic fields. It turns 
out that the unique coexistence of weak ferromagnetism 
and chiral modulations enables the occurrence of new 
types of incommensurate structures and specific local- 
ized structures including so called magnetic vortices or 
skyrmions which are generally unstable in other classes 
of magnetic materials. 



II. THE MODEL 



A. Phenomenological energy 



The tetragonal antiferromagnet Ba2CuGe20vli3'lI_l 
(space group P42ito) belongs to the crystallographic 
class D 2 d and K2V3O8LJ (space group PAbm) to Ca v . 
The magnetic (free) energy within a continuum descrip- 
tion consistent with the symmetry and the two-sublatticc 
magnetic structure of these antiferromagnets can be de- 
rived by the standard approach to phenomenological 
theory.El At temperatures sufficiently below the ordering 
temperature the vectors of sublattice magnetization M.; 
(i = 1,2) do not change their modulus. In this prac- 
tically important case defined by neglecting the para- 
process, the vectors Mi have only orientational degrees 
of freedom and can be described by the unity vectors 
mi = Mi/M s , where M s = |Mi| = |M 2 | is the sublattice 
saturation magnetization. For tetragonal antiferromag- 
netic crystals, the two-sublattice model described by the 
unity vectors mi and m 2 yields the magnetic energy in 



the following form: 




-h • (mi + m 2 ) 



> J / 2 , 2 \ a' 

- 2 ( m iz + m 2z) ~ P m lz m 2z 
-d(mi x m 2y - m 2 xmiy) + w D ^j dV. (3) 

This includes inhomogeneous (a, a') and homogeneous 
(A) parts of the exchange coupling and the interac- 
tion energy with the external field h. The next two 
terms describe uniaxial second-order magnetocrystalline 
anisotropy with constants (3, /?', where the z-axis is taken 
along the tetragonal axis of the antiferromagnets. The 
homogeneous part of the Dzyaloshinskii-Moriya interac- 
tion with constant d is responsible for weak ferromag- 
netism with small magnetic moments in the basal plane. 
Finally the energy contribution wd includes Lifshitz in- 
variants of type (g). The functional form of wrj depends 
on the crystal symmetry and will be specified later. 

The next terms in a systematic expansion of the energy 
for a two-sublattice antiferromagnet are much weaker 
fourth-order terms of the magnetocrystalline anisotropy, 
including uniaxial parts with terms mf z , m 2z , vn\ z rn\ zl 
and a magnetic anisotropy in the basal plane (XOY 
plane) composed of x- and y-components of the vectors 
m t . The former are important in close vicinity to some 
reorientation transitions and the latter is responsible for 
small variations of magnetic structures when the mag- 
netic field is rotated in the basal plane. These secondary 
effects are omitted in this contribution dedicated to the 
principal magnetic properties of the system. We also 
neglect the stray field contribution in the total energy 
(H) because, due to the antiparallel alignment of mag- 
netic moments in antiferromagnets, stray fields are much 
weaker than in ferromagnetic crystals. They, however, 
play a crucial role in stabilizing multidomain structures 
in the vicinity of field-induced reorientation transitions.113 

The functional (^) includes all leading interactions 
in an uniaxial two-sublattice antiferromagnetic crystal. 
Here, we briefly list several special cases of the model (||) 
which describe important special classes of antiferromag- 
netic systems. 

I. ) d = 0, wrj = 0. Collinear antiferromagnets . The vast 
group of these antiferromagnetic materials includes sueh 
well-studied species as CuC1 2 -2H 2 OEI, MnF 2 B, Cr 2 3 E3, 
GdA103t3 (see for furthet-cefereiices and review of their 
magnetic properties Refs.E3 andcJ). 

II. ) d ^ 0, wd = 0. Antiferromagnets with weak fer- 
romagnetism. In this case the energy (|3|) describes anti- 
ferromagnetic crystals with homogeneous Dzyaloshinskii- 
Moriya interactions resulting in weak ferromagnetism 
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or collinear antiferromagnets with "hidden" weak ferro- 
magnetism. Among many others, this, group includes 
MnCO3,0 orthoferrites]23 manganitesj2a and the most 
popular and well-studied weal^fe^mapiaffnetic antiferro- 
magnet hematite, i.e. a-Fe2O3.tr0OE3'E2l 

III. ) d — 0, wd 7^ 0. Chiral helimagnets. This case is 
realized in the cubic helimagnets discussed ahpije and in 
other noncentrosymmetric magnetic systems.LrO Usually 
interactions of type (phpSiabilize modulated chiral struc- 
tures in these materialsllj. 

IV. ) d 0, wd 7^ 0. Chiral helimagnets with weak fer- 
romagnetism. The model (||) represents previously un- 
known systems where both, homogeneous and inhomo- 
geneous, Dzyaloshinskii-Moriya interactions are opera- 
tional and not forbidden by any other additional sym- 
metries. As will be shown in this paper, this unique co- 
existence of a mechanism for weak ferromagnetism and 
chiral coupling leads to specific modulated states with 
magnetization oscillations. Thus, in the following we will 
generally consider systems d ^ and Wd 7^ 0. 



B. Simplified model and basic equations 

It is convenient to use linear combinations of the sub- 
lattice magnetization m^, namely, the vector of total 
magnetization m = (mi + m.2)/2 and the staggered mag- 
netization (or vector of antiferromagnetic order) 1 = 
(mi — rri2)/2 as internal variables of the system. Be- 
cause of I nii| = 1 these vectors satisfy the constraints 
m • 1 = 0, m 2 + l 2 = 1. 

In most antiferromagnetic crystals the exchange cou- 
pling is much stronger than other internal interactions. 
Strong magnetic fields of order A destroy antiferromag- 
netic order and orientate the sublattice magnetizations 
m^ parallel to each other (a so called spin-flip transition 
into the "paramagnetic" phase with |m| = 1, 1 = 0). For 
most investigated antiferromagnetic systems these "ex- 
change" fields are extremely large. Practically attainable 
values of magnetic fields usually only slightly distort the 
antiparallel arrangement inducing states with the total 
magnetization much smaller than unity. The hierarchy 
for the strength of interactions, A 3> d,/3,(3', and the re- 
lations for the internal parameters, m«l and per- 
mit to considerably simplify the energy (|J) by excluding 
gradients of m and takingiato account only the following 
terms (for details, see e.gra) 



quadratic in the components of 1 have the following form 



W =J{ A 



3 

E 

1=1 



75— + Am 2 

OXi 



2m- h 



-2d (m x l y 



Tflylx) 



Bl 2 z +w D \dV, (4) 



where A = a — a' and B — (3 — (3' . Functional forms of 
wrj for all noncentrosymmetric crystallographic classes 
have been derived in Ref. |3^. In particular, for the anti- 
ferromagnets under consideration the Lifshitz invariants 
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dL 
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class D 2 d ■ wd — D l z — l x -^~ + h^r 1 - - lyTr- 



dy dy 



dx 



dx 



(5) 



class C nv : wd = D [ I 



dl x 
' dx 



ox ay Oy 



(6) 



The homogeneous part of the Dzyaloshinskii-Moriya in- 
teraction in (jij) includes in-plane components of m and 
1. They originate from the z-component of the vector 
product ([j]). Writing the Dzyaloshinskii vector in ([|) 
as a sum of two parts proportional to D and d, this 
contribution to (0) can be derived from the vector di- 
rected along the tetragonal axis d = (0,0, d). On the 
other hand, the Lifshitz invariants (j^), (|J) can be de- 
rived by an expansion for the in-plane components of the 
vector product ([I]) considering the contribution due to 
the vector D = (D,D,0). The terminology in this field 
is not yet fixedly formulated. Here, following Ref. |3l], we 
will call energy contributions given by Lifshitz invariants 
Dzyaloshinskii or chiral interactions to distinguish them 
from the homogeneous part of the Dzyaloshinskii-Moriya 
interaction. 

Independent minimization of the energy (^) with re- 
spect to m leads to the following result 



m = — [n x (d + nx h)]/A , 



(7) 



where n = 1/|1| is the unity vector parallel to the stag- 
gered magnetization vector. After substitution of (fjj) into 
the energy (Jij) one obtains the energy to leading approx- 
imation as a function of the vector n 



di 



A 



[{h x + dn y f + (h y - dn x f - (h • n) 2 ] 



-Bn 2 z + w D (n) } dV . 



(8) 



where wo(n) is determined by (|) or (^). The energy (||) 
and (0) were derived from ([|) by ignoring the parapro- 
cess and assuming weak total magnetization, |m| <C 1, 
implying |1| ~ 1. Both assumptions are fulfilled in most 
realistic cases of interest. Thus, the energy (||) describ- 
ing the orientation of the staggered magnetization can be 
considered as general phenomenological description for 
realistic uniaxial two-sublattice antiferromagnets. The 
functional (0) is related to so-called nonlinear cr-models 
which are basic subjects in the theory of solitons and 
which arc intensively studied in mathematical and theo- 
retical physicsEj 

The energy contributions (^) and (JsJ) can result in 
states with modulations in the basal tetragonal plane, i.e. 
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with a propagation vector in the XOY-plane. Its actual 
direction is selected by some small in-plane anisotropy 
contributions that are neglected here (see discussion on 
expansion (|J) above). On the other hand, there are no 
interactions in our systems violating homogeneity along 
the tetragonal z-axis. Hence, we infer that within the 
model (||) the most general solutions are inhomogeneous 
only in the basal plane but homogeneous along the z- 
axis. It is convenient to write the vector n(x, y) and the 
magnetic field h in spherical coordinates: 



n = (sin # cos ^, sin sin i/;, cos 6 1 ), 

h = (h sin ( cos 77, h sin £ sin 77, h cos £) . 

In these variables the total energy is given by 
W =L Z J{ (^f+sin^g 



(9) 




+u>d + w > dx dy . 



(10) 



where the integration with respect to z was performed 
for a system with linear size L z , x± = x, xi = y. The 
Lifshitz invariants are given by 
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for D 2 d ■ 



w D 
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cos ip 
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(12) 



sin cos ( 



, Otp , 8i\) 

sin^- cosi/i— — 

Ox Oy 



for C n 



and the energy term w does not depend on spatial deriva- 
tives: 

A* = —(\B - d 2 — h 2 cos 2 C) cos 2 6 

-(h 2 sin 2 C + d 2 ) + h 2 sin 2 ( cos 2 (ip - 77) sin 2 9 
—2dh sin £ sin 9 suv(ip — 77) 

+h 2 sin £ cos £ sin 29 cos(ip — rf) . (13) 

The functional (|l(]) with (|l^) provides the basic expres- 
sion for the total energy of uniaxial two-sublattice antifer- 
romagnets belonging to crystallographic classes without 
inversion symmetry. By inserting the appropriate Lif- 
shitz invar iant from (|TT|) or (O) for wd, the functional 
describes the magnetic energy of the two tetragonal crys- 
tals of interest here. 



III. THE PHASE DIAGRAM OF EQUILIBRIUM 
SOLUTIONS 



The equilibrium distributions of 9(x, y) and ip(x, y) are 
determined by solving a set of equations minimizing the 



energy (]l(]) . Depending on the values of the phenomeno- 
logical constants in the energy ( |l~0| ) and the components 
of magnetic fields different spatially homogeneous and 
modulated phases can be realized in the system. Due to 
isotropy of the model in the basal plane only the com- 
ponent of magnetic field along the tetragonal axis (h z ) 
and the value of its projection onto the basal plane h± 
are of importance. A reduction of the number of control 
parameters is obtained by rescaling the spatial variables 
and the energy. We use the following units for lengths, 
magnetic field, and the strength D of the inhomogeneous 
Dzyaloshinskii-Moriya interaction 



x = VAX/\K\, 



D = -y/A\K\/X, 




(14) 



introducing an effective anisotropy constant K acting 
on the staggered magnetization that is comparable to a 
constant of uniaxial anisotropy in ferromagnets. In cen- 
trosymmetric antiferromagnets (D = 0) and zero exter- 
nal field, the collinear state with staggered magnetization 
1 along the tetragonal axis and m = is the ground state 
for K > (easy- axis system); for K < the vectors 
1 and a nonzero m lie in the in the basal plane (weak 
ferromagnetic states). The characteristic length xo is of 
the order of the effective size of an isolated domain wall 
between the homogeneous states at zero field. In uniax- 
ial ferromagnetic materials the corresponding expression 
for an aHtrinsic length is known as exchanye or Block 
length. c3 The characteristic field ho is the so-called spin- 
flop field. Finally, the parameter Do equals the lowest 
value of the (Dzyaloshinskii) constant D that stabilizes 
modulated states at zero field (see below). In these re- 
duced units ( |i~4| ) , the energy ([l0]) includes as independent 
parameters the rescaled constants K involving d and D 
as well as the two components of the applied field (h z , 
h±). Thus, these parameters span a four-dimensional 
phase-space for the solutions. 

Before giving the detailed analysis, let us point out 
some general features of the possible magnetic configura- 
tions in this system. The equilibrium magnetic structures 
are governed by two opposing tendencies. The rotation 
of the staggered vector 1 with propagation vectors in the 
basal plane and an appropriate sense of rotation leads to 
negative values of the invariants ( p4[ ) , (|l^) . An unlimited 
reduction of the pitch for this winding of the staggered 
magnetizations would lead to infinitely negative values 
of this Dzyaloshinskii energy. This is counter-acted by 
the inhomogeneous part of the exchange energy in (|10| ) 
providing the "stiffness" of the magnetic structure. In 
isotropic systems, i.e. w=0 for expression (|l3|), the ratio 
of these competing energy contributions yields the optL 
mal period for the spiral, which is of the order of A/Dxi 
Such chiral modulations with uniform rotation are ott. 
served in low-anisotropy systems as cubic helimagnctsLJ 
or in hexagonal chiral magpets with in-plane rotation of 
the magnetization vectors.t3 The uniform rotation of 1 
in spirals is disturbed by anisotropic interactions and/or 
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by application of a magnetic field, i.e. the energy terms 
included in ([T^). These interaction terms result in pre- 
ferred directions for the staggered vector 1 corresponding 
to the minima of the energy density (|l3|). Hence, they 
distort the chiral modulations and may even suppress 
them by forcing the staggered magnetization to point 
fixedly into "easy" directions. Thus, chiral modulations 
may occur only beyond a certain threshold: the interac- 
tions (||) or (||) must be strong enough to overcome the 
anisotropic energy contributions suppressing modulated 
states. Below this threshold the system takes on the ho- 
mogeneous states which are determined by minimization 
of the energy ( [l0| ) with D = A = 0. In the following sub- 
sections we will demonstrate this competition between 
homogeneous and inhomogeneous states for our model in 
detail. 

Another important property of the system is related 
to the role played by the homogeneous Dzyaloshinskii- 
Moriya interaction: in-plane components of the staggered 
vector induce corresponding non-zero components of the 
total magnetization vector m (weak-ferromagnetic mo- 
ments) according to eq. (|t]). Thus, collinear antiferro- 
magnetic states, m = 0, exist only when the staggered 
vector 1 is parallel to the tetragonal axis. All magnetic 
structures with a vector 1 deviating from this direction 
perforce have a locally nonzero magnetization m. In zero 
field, these magnetic moments are in the basal plane. 
Helix-states of 1 are accompanied by magnetization com- 
ponents oscillating in sign and with the same period as 
the antiferromagnetic modulations. This peculiar mech- 
anism leading to modulated antiferromagnetism and a 
related magnetization m may become operational and 
important, even if the ground-state is not modulated, for 
the more general case h ^ (see below). 

In the general case the four dimensional phase dia- 
gram (d, D, h z ,h±) includes regions with different mod- 
ulated and homogeneous states separated by "hypersur- 
faces" corresponding to different phase transitions. As 
remarked above, there are numerous experimental and 
theoretical results on magnetic properties of centrosym- 
metric antiferromagnets (model.j£2l with wd = 0) de- 
scribing weak ferromagnetism.Ejt_E_l The phase space of 
the control parameters in this case (d, h z ,h±) was found 
to have a very complex topology and, depending on the 
orientations and the relative strengths of the vectors h 
and d, ajiumber of non-trivial transitions occur in these 
systemsBS This phase diagram can be considered as a 
"cross-section" given by the three-dimensional "hyper- 
plane" D = through the general (d, D, h z , /ijj-space 
investigated here. On the other side in chiral antiferro- 
magnets described by the case d = for our model, both 
the spin arrangements in a spiral and the correspond- 
ing propagation directions are found to be very sensitive 
to the orientation and strength of the applied fieldEJ In 
the general case of nonzero values of the constants d and 
D, there is an even wider variety of homogeneous and 
inhomogeneous solutions characterized by complex non- 
collinear magnetic structures and variable directions of 



propagation vectors. Clearly, the full set of the homoge- 
neous and inhomogeneous states corresponding to the en- 
ergy ( |l0| ) together with (|ll]) to ([l3]) are of general interest 
and well-worth of further investigations. Here, we restrict 
ourselves to the approximation valid for antiferromag- 
netic systems with small magnetic anisotropy and their 
specific hierarchy of the interactions. The physically ex- 
pected relation d^> B turned out to be valid in all known 
systems with weak ferromagnetism and is based on the 
common rclativistic origin of both magnetic energy con- 
tributions in these (low-anisotropy) systems. It is specif- 
ically valid for systems in which magnetism is due to 
d-electrons. In these systems, the Dzyaloshinskii-Moriya 
interactions generally overcome the magnetocrystalline 
anisotropy. This is expected to apply also in noncen- 
trosymmtric antiferromagnets with d-electron magnetism 
for which the values for d are still unknown. This allows 
to narrow considerably the range of physically meaning- 
ful control parameters in our model. 

As already discussed above, the vector d induces a 
total magnetization in the basal plane which tends to 
orientate itself parallel to the in-plane components of an 
applied field h± . Correspondingly, the staggered magne- 
tization is rotated into the plane perpendicular the h± 
direction. The stronger the values of d, h± and of the 
in-plane components of 1 the stronger is this effect. We 
may assume that the deviation of the staggered magneti- 
zation from this plane perpendicular to the h± direction 
e = \ip — r) — 7r/2| is small. By optimizing the energy 
(|l3| ) with respect to e one obtains 



| V - 1? - tt/2| 



h cos £ cos 9 



d + h sin C sin t 



(15) 



providing the consistency criterion for our assumption 
e « 1. This approximation is valid always for d 3> h 
and generally in a broad range of the orientations for 
the vectors h and 1. This includes almost all physically 
interesting cases. In the following analysis, we assume as 
central approximation that the staggered magnetization 
is always restricted to the plane perpendicular to i.e. 
e = 0. By substituting ip — r\ = tt/2 the energy density 
( |l3| ) can be simplified and reduced to the following form 
using the scaled quantities (14) 



$(0) with 



\K\ 



$(0) = sgnJf ( 1 - — cos^ C ) (sinf 



\ 2 



dh sin £ 



K — h 2 cos 2 £ 



(16) 



(Here, we drop constant terms in w, i.e. those indepen- 
dent of 0.) In the following subsections, we investigate 
spatially homogeneous phases, helical phases, and their 
respective stability limits with the approximation (|l6|). 
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FIG. 1: Basic spin configurations described by the staggered magnetization 1 and the total magnetization m (for clarity, 
|m| -C |1| — 1 is not obeyed in the pictures) in homogeneous states of easy-axis systems (K > 0). 

(a) and (b) external field h along the tetragonal axis, (a) Collinear or antiferromagnetic (AF) phase with 1 || OZ and zero 
magnetization exists in the applied field along the tetragonal axis < h < ho- 

(b) Spin-flop phase with m || OZ and with 1 in the basal plane is an important particular case of the weak ferromagnetic (WF) 
phase for h > ho along the tetragonal axis. 

(c) and (d) magnetic field with oblique direction. The low symmetry canted phase (c) exists for v < 1 (jl^) and continuously 
transforms into the WF phase (d) at the critical line v = 1. In the canted phase the staggered magnetization is in the plane 
perpendicular to the the in-plane component of the applied field and all components of the vector m have generally nonzero 
values (cf. eg. (0). (d) General case of the WF phase. The staggered magnetization 1 lies in the basal plane, the component 
of the total magnetization along the tetragonal axis (m z ) is induced by the corresponding component of the applied field, the 
in-plane components are due to the homogeneous Dzyaloshinskii-Moriya interaction and in-plane components of magnetic field. 
In easy-plane systems (K < 0) only the WF phase is stable for all values of the applied field. 



A. Homogeneous states 

The homogeneous states are described by the be- 
haviour of the energy functional (fl6|). Depending on the 
sign of K the energy (|l6| ) describes two different types of 
antiferromagnetic ordering. 



1. K > Easy-axis system 

At zero field and in a magnetic field along the tetrag- 
onal axis for h < ho the antiferromagnetic phase with 
1 || z-direction and m = has the lowest energy. This 
magnetic structure is sketched in Fig. |l|(a). At the field 
ho = \[K the vector 1 "flops" down onto the basal plane. 
This is a so-called spin- flop transition. In the resulting 
spin-flop phase with 9 = tt/2, the total magnetization un- 
der influence of the homogeneous Dzyaloshinskii-Moriya 
interaction is slightly inclined from the tetragonal axis 
(Fig. |l|(b)). In the region where |m| <C 1 the compo- 
nents of m, m z — h/X, \m±\ — d/X, are obtained from 



(Q). In the region h$ < h < X the total magnetiza- 
tion increases linearly with increasing field and finally at 
the "exchange" field h e x — X the spin-flop phase con- 
tinuously transforms into the saturated "paramagnetic" 
phase with |m| = 1, 1 = by a spin-flip transition. 
Note, in the spin-flop phase the magnetic state has an 
infinite degeneracy with respect to rotation of the vec- 
tor 1 around the tetragonal axis. In-plane anisotropy 
reduces the degeneracy to certain preferable directions 
related by symmetry in the basal plane. E.g., in the 
case of fourth-order tetragonal anisotropy there are two 
mutually perpendicular directions of "easy" magnetiza- 
tion. In an increasing magnetic field deviating from the 
tetragonal axis the staggered magnetization 1 rotates to 
the basal plane in the plane perpendicular to the pro- 
jection of h onto the basal plane. The angle 9 between 
the vector 1 and the z-axis is sin (9 = v for v < 1 Jl6|). 
We name this state with a finite angle between staggered 
magnetization and z-axis canted phase (Fig. [l](c)). Fi- 
nally, at the critical line h c (h±, h v ) (Fig. ||) where v from 
expression ( |l6| ) attains the critical value v = 1, a phase 
transition occurs into a phase with the staggered vector 



weak ferromagnetic phase 



h c (h\,hy) 




FIG. 2: (hj_,hz) phase diagram for the homogeneous states 
in easy-axis systems [K > 0). The low symmetry canted 
phase transforms into the WF phase by a second order phase 
transition at the critical line h c (h±,h z ) determined by the 
condition v = 1 in (|l^). Scales are given by ho = K 1 ^ 2 and 
h° x = K/d. 



lying in the basal plane, sin 9 = 1, and perpendicular to 
the applied field (Fig. |l|(d)). This is a weak ferromag- 
netic (WF) phase. The phase-diagram for this transition 
between canted and WF phase is depicted in Fig. g. The 
total magnetization is deduced by substituting the equi- 
librium values of n — 1/|1| into eq. (Q). Assuming that 
the applied field is in the XOZ-plane, then the staggered 
magnetization rotates in the KOZ-plane. From (uh the 
following expressions for the magnetization components 
result 

m x = (/isin£ + <isinr9)/A, 

m y = —hcos(s'm9cos9/\, 

m z = ft, cos £ sin 2 9/X , (17) 

where 9 = arcsini/ for the canted phase and 9 = tt/2 in 
the weak-ferromagnetic phase. 

The relative orientation of the vectors 1 and m is fixed 
by the sign of the constant d. This handedness of the 
magnetic structures reflects the chiral character of the 
interaction (Q) and leads to the nonequivalence of ener- 
gies for the states with antiparallel values of 1 in oblique 
magnetic fields. One can understand this considering the 
shapes of the potential <I> (|l6| ) with applied fields (Fig. ||). 
In the general case (Fig. 131(a)) of an oblique field, i.e. an 
applied field deviating from the tetragonal axis, the anti- 
ferromagnetic phase with 9 = %n transforms into the WF 
phase (h > h c ) via the canted phase. From the potential 
profile for this canted phase (Fig. ||(a), < h < h c ) one 
immediately sees for a given state that the correspond- 
ing state with antiparallel orientation of 1 has a different 
energy and generally is not an equilibrium state. Thus, 
the symmetry between states with antiparallel 1 peculiar 
to ordinary antiferromagnetic phases is violated. The 
stable states in this canted phase are separated by two 
types of potential barriers. At the transition into the 
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FIG. 3: Schematic evolution of the potential profiles <&(#) un- 
der influence of the applied field for easy-axis systems (K > 0) 
in an oblique field (a) and in a field in direction of the tetrag- 
onal axis (b). 

In (a) for the canted phase two different barriers occur be- 
tween equivalent equilibrium states. Different domain walls 
correspond to these barriers: DW I between state (1) and (2), 
and DW II between (1) and (3). 



WF phase the lower potential barrier disappears and the 
higher separates states with 9 = ir/2+2mr. In a magnetic 
field along the tetragonal axis (Fig. ||(b)) the potential 
barriers in $ ( Jl6| ) disappear as the field approaches the 
spin-flop field from both sides, i.e. in the antiferromag- 
netic {h < ho) and in the spin-flop phases (h > ho). This 
means that at the transition field h — hg the potential 
barriers between coexisting antiferromagnetic and spin- 
flop states are anomalously low and are determined by 
the values of the fourth-order anisotropy. 



2. K < 



Easy-plane system 



In the ground state for easy-plane systems K < the 
staggered magnetization lies in the basal plane with a 
spontaneous magnetization, |m| = d/X, perpendicular to 
the vector n. Therefore, this is a weak-ferromagnetic 
phase. The behaviour of the system under influence of 
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(a) class D 2d 




FIG. 4: Basic modulated structures (a) - helicoid for systems 
with D2d symmetry; (b) - cycloids for antiferromagnets with 
Cnv symmetry. 
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FIG. 5: Typical phase portrait of the solutions for the equa- 
tion (119) in the canted phase. The separatrix curves between 
open (continuous lines) and closed (dashed lines) orbits are 
highlighted by a thick line. 



a magnetic field is similar to that in the easy-axis sys- 
tem (section [II A 1) for h z larger than the spin- flop field, 
namely the vector 1 is oriented perpendicular to h, and 
the total magnetization gradually increases in increasing 
field. 



B. Helical structures 

The equations minimizing the functional ([lO]) also per- 
mit solutions with chiral modulations propagating in 
the basal plane. First we consider structures modu- 
lated along a certain fixed direction in the basal plane 
and homogeneous perpendicular to this direction. This 
yields one-dimensional spirally modulated states com- 
prising helicoids and cycloids. In the absence of in- 
plane anisotropy all propagation directions are equiva- 
lent. The structure of these modulated states depends on 
the crystal symmetry which manifests itself in different 
functional forms of the Lifshitz invariants (see (jjj), 
For antiferromagnets belonging to the crystallographic 
class £>2d the staggered vector 1 rotates in the plane per- 
pendicular to the propagation direction, i.e. ip = ixjl. 
These states are helicoids (Fig. ||(a)). This rotation of 
1 reminds the behaviour of the magnetization vector for 
Bloch walls in ferromagnets (seet3). In the case of C nv 
symmetry, 1 rotates in the plane formed by the tetrago- 
nal axis and the propagation direction (-0 = 0) forming 
cycloids (Fig. f|(b)). This is akin to Neel domain walls 
in ferromagnets. The rotation in the spirals has a fixed 
sense determined by the condition that the inhomoge- 
neous Dzyaloshinskii-Moriya energy must be negative. In 
both cases, the spirals are accompanied by oscillations of 
the magnetization m perpendicular to the plane of rota- 
tion according to ^ (Fig. |]). In zero field, m = d sin 6>/A. 
Under the influence of an applied magnetic field the spi- 
rals orientate in such a way that the rotation of 1 oc- 
curs in the plane perpendicular to the projection of h 



onto the basal plane. As above for the homogeneous 
structures, we assume that the magnetic field lies in the 
XOZ-pl&ne and the vector n in the YOZ-pl&ne. For 
helicoids (r>2d symmetry) the propagation direction is 
along the x-axis and for cycloids (C nv symmetry) along 
the y-axis. The spatial coordinate along the propagation 
direction is measured in reduced units of xq according to 
({Til), £ = Xi/xo. In these reduced units the energy func- 
tional © for one-dimensional modulations assumes the 
form 



W = 




AD 



tvDq \d£_ 



(18) 



where the multiplicative constants due to the integrations 
in directions of z and perpendicular to £ are absorbed in 

W. $(6>) is given by (|l6|) and D Q by @. The first 
integral of the Euler equation for the functional (|l8|) is 
readily derived, 



- $(6») =E, 



(19) 



where E is an integration constant. In passing, we re- 
mark that the Euler equation with the potential $(#), 
from ( |i"6| ) is related to the double sine-Gordon equation.EJ 
Typical phase trajectories of (|l^) in the (#£, 6 i )-phase 
plane are plotted in Fig. |E| (here, we use the abbrevia- 
tion = dd/d£ ). The separatrices obtained for E = 
cross each other in the points corresponding to the min- 
ima of the function $(#). They divide the phase plane 
into regions with closed (E < 0) and open (E > 0) 
trajectories (Fig. ||). The closed trajectories correspond 
to alternating rotation of the staggered vector and obvi- 
ously are not of interest for our model as they describe 
inhomogeneous states with alternating sense of rotation 
that do not minimize the energy related to the inhomo- 
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FIG. 6: Series of solutions of ( |fg| ) #(£, E) corresponding to 
the open orbits in Fig. 5. Solutions shown by thinner lines 
correspond to open orbits with E > close to the separatrix 
curve of Fig. 5 (not shown there). 



geneous Dzyaloshinskii-Moriya interactions. The mod- 
ulated states with fixed rotation sense are described by 
open trajectories. The integration of eq. ( |l9|) with $(#) 
from ( |l6| ) yields the set of solutions #(£, E) parametrized 
by the constant E (Fig. ||). These solutions can be ex- 
pressed analyticallyr-as certain cumbersome combinations 
of elliptic functions^ (see, e.g., Refs. |^, f30p4| ) . Here, for 
simplicity, we derive representative solutions by direct 
numerical integration of (|l9|). Using ( |l9| ) the energy den- 
sity w averaged over a period S of a modulated state can 
be written as functions of the parameter E 



1 



[E + 2$(6)]d6 2ttD 



S(J5) ,o 

271 



5(25) = / 
Jo 



d6 



E 



3(E)- 



(20) 



Note, eq. ( ]l9| ) does not involve contributions from the 
inhomogeneous chiral Dzyaloshinskii-Moriya interactions 
and, thus, its solutions (Fig. ^) do not depend on these 
chiral interactions. Therefore, these solutions (Fig. ||) 
have the same functional form as those for the corre- 
sponding model with D = (centrosymmetric systems). 
However, the energy w (|20| ) of the system depends on the 
contribution from the Lifshitz invariants. This energy has 
different values for different integral curves 6(£,E). The 
equation dw/{dE) = to derive the optimal values E 
can be reduced to the following form 



/ d6^$(6) + E = 4D/D . 
Jo 



(21) 



Hence, the spiral structure described by the integral 
curve 0(£, E) obeying (|2l]) corresponds to the equilib- 
rium magnetic structure realized in a noncentrosymmet- 
ric system, where inhomogeneous Dzyaloshinskii interac- 
tions wd are operational. From the solutions of eqs. jl9| ) 
and ( pf| ) the other equilibrium parameters of the spiral 



structures are readily calculated. In particular, eq. (|2£ 
yields the period of the structure 2, and the oscillating 
components of the vector m are expressed via as func- 
tions of #(£, E). Depending on the ratio D/Dq for the 
relative strength of the chiral interactions the modulated 
structures display the following characteristic evolution: 
For strong chiral interactions D/Dq 3> I the influence 
of the energy (|l6| ) is negligible and in the equilibrium 
states the staggered vector rotates with an essentially 
fixed "velocity" 0^ corresponding to the phase trajecto- 
ries with = 2D/(ttD ). Perturbations of the uniform 
rotation for the spirals are related to the shape of the 
potential profiles (cf. Fig. [I| and the open trajectories 
in Fig. ||). Thus, the functional dependencies of 0(£,E) 
contain important information on internal interactions 
of a system. When values of the parameter are smaller, 
D/Dq ~ 1 the influence of the potential $(0), which de- 
termines preferable orientations in the crystal, violates 
uniform rotation of the staggered magnetization 1 in a 
spiral. Further weakening of the chiral interactions leads 
to "pinning" of 1 along certain ("easy") directions and 
squeezes the regions with "disadvantageous" orientations 
of the vector 1. This tendency results in the formation of 
structures consisting of large domains with homogeneous 
states separated by narrow transition regions in which 
the vector 1 rotates from one easy direction to another 
similar to domain walls. Finally, for a certain critical 
value of D the modulated phase is transformed into the 
homogeneous state. This transition is signalled by an un- 
limited growth of the period for the modulated state. In 
the phase space, states at this transition into the homo- 
geneous state correspond to the separatrix (Fig. ||) which 
describes a set of isolated domain walls, i.e. walls with 
infinite separation between them. The finite stiffness of 
the exchange interaction prevents a complete annihila- 
tion of the domain walls and they may exist with finite 
thickness within homogeneous states as metastable topo- 
logically stable objects acting as nucleation centers dur- 
ing a reversal transition from the homogeneous into the 
modulated state. 

Finally, it should be stressed that eqs. ([jl]) and 
( ^l| ) provide general and rigorous solutions for one- 
dimensional modulated structures in magnets with Lif- 
shitz invariants of type @ with arbitrary functional form 
for the potential < &(6') in the functional fll8|). The above 
described evolution of the modulated states is not re- 
stricted to any particular form of <&(#); the qualitative 
picture of this evolution rather has universal character 
for physically reasonable choices for $(0). Particular 
cases for such chiraLspirals have been investigated start- 
ing from the paperu for.-s evfflra l groups of helical ferro- 
and antiferromagnets.y'EJ'L^EZI 



C. Stability limits of the modulated states 

At the transition into homogeneous states the chiral 
spirals disintegrate into a system of noninteracting pla- 
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nar domain walls. Such a transition can be found by 
comparing the energy of the chiral spiral to the energy 
of domain walls separating regions with different homo- 
geneous states of the system.HEll Below in section [f"v| , 
we will discuss domain walls as topological defects in 
the magnetic system which require an excitation energy. 
Here, we are concerned with the competition between 
homogeneous and modulated equilibrium states. Then, 
a gain of energy through proliferation of domain walls in- 
dicates the instability of homogeneous states compared 
to a modulated state. 

Let us consider a planar isolated domain wall be- 
tween two infinitely extended regions with different spa- 
tially homogeneous magnetic structures that are de- 
scribed by the functional (|l6|), i.e. equilibrium states 
in Fig. pi The equilibrium structure of this isolated 
wall is determined by solving (091) with the bound- 
ary conditions 9 |^ = ±oo= 0i,2> \d9/d(,)^ = ±oo = 0, 
where 9\ t 2 are homogeneous configurations determined 
by $oo = min[$(0)] = $(0i) = ^(62) for the functional 
(jig). The direct integration of (^9|) yields the following 
results for the dependence|-of-Q(^) in the wall and for the 
domain wall energy a (seeM) 



d0' 



8l vw) - $c 



7T „ 

2 D ° 



A/[$(0)-$oo]d0 ± £>|0i - e 2 



(22) 



(23) 



K<0 

weak ferromagnetic phase 
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spin-flop phase 
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The function [$(6) — $00] is the deviation of the energy 
density @ fr om the minimal value Woo corresponding 
to the homogeneous states in adjacent domains. The 
first term in ( |23] ) is positive and represents increased en- 
ergy contributions compared to those of the homogeneous 
states 9\, #2- This increase is due to inhomogeneous 
exchange interactions and interactions included into the 
functional $ ([if])- These defect energies are typical for 
the energy of magnetic domain walls £3 The second term 
is specific for noncentro symmetric systems. Its sign is de- 
termined by the rotation sense of the staggered vector n 
in the domain wall. Clearly, for any sign of the constant 
D there exists a rotation sense of the staggered magne- 
tization leading to negative values for this energy con- 
tribution and, consequently, to a decrease of the domain 
wall energy. For sufficiently strong inhomogeneous chiral 
Dzyaloshinskii-Moriya interactions the total energy of the 
domain wall may be negative compared to the energy of 
homogeneous states. This manifests an instability of the 
homogeneous state with respect to chiral modulations. 
As already discussed above, such a transition takes place 
for coefficients D larger than a certain threshold value 
necessary to overcome the positive energy contribution 
for inhomogeneous states due to the conventional mag- 
netic interactions. 

The wall energy (^3|) can be expressed via the height of 
the potential barrier, A<f> = max[<&(0)] — $oo, that sepa- 



FIG. 7: Magnetic phases in the magnetic field along the 
tetragonal axis in dependence on strength D of the inhomo- 
geneous Dzyaloshinskii-Moriya interactions (a) K < 0; (b) 
K > 0. Note that for easy-axis antiferromagnets (b) in the 
vicinity of spin-flop field the modulated phases exists at arbi- 
trarily small values of D. 



rates the equilibrium states 9\ and #2 ■ This can be writ- 
ten in the following form 



(24) 



where 7 is a numerical factor determined by the average 
value of the integrand in (|23|). These results have clear 
physical meaning. The higher the energy barrier A<i> the 
stronger the chiral interaction necessary to overcome it 
and to stabilize modulated states. As was shown above, 
the potential profile $(0) ( |l6| ) strongly depends on the 
strength and direction of the applied field (Fig. ||) . Cor- 
respondingly, the critical values D c of the Dzyaloshinskii- 
Moriya constant D for transitions between homogeneous 
(D < D c ) and modulated chiral states (D > D c ) vary 
strongly with an applied magnetic field. The critical sur- 
faces for these transitions in parameter space are given 
by the equation er = using the functional (nw. For 
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easy-plane systems (K < 0) the critical surface is 



D 



hi 

-J - 1 

K 



(25) 



This equation also describes the boundaries of the mod- 
ulated states for easy-axis systems (K > 0) in magnetic 
fields larger than the spin-flop field (h z > ho ). For lower 
fields (h z < ho ) the equation 




V 7 ! -v 2 



v arcsm v 



(26) 



gives the transition into the canted phase (y < 1). Fi- 
nally, the equation 




v arcsm 



(27) 



describes the transition into the weak-ferromagnetic 
phase (y > 1). In particular, at zero field the critical 
value D c equals Do- Thus, this constant is the lowest 
value for the Dzyaloshinskii constant D to induce mod- 
ulated ground states. The (D,h z ) phase diagrams for an 
applied field in direction of the tetragonal axis are shown 

in Fig. [?] for the two cases K > . 

For easy-plane systems (K < 0) the critical surface 
D c (h±, hz) has a minimum in the origin with -D c (0, 0) = 
Dq and monotonically increases with increasing magnetic 
field for any direction (Fig. ||(a)). When D < D the chi- 
ral interactions are too weak to overcome the pinning due 
to the uniaxial easy-plane anisotropy. Then, the system 
exists in the homogeneous state with the staggered mag- 
netization in the basal plane and a weak spontaneous 
magnetization (WF phase). For D > Do the WF phase 
becomes unstable. Under the influence of the inhomoge- 
neous chiral interactions the vector n "escapes" from the 
basal plane and a chiral helix is formed. We add a remark 
about the peculiarity of this type of helix. In known easy- 
plane systems with helical structures the magnetization 
(or the staggered magnetization in the case of antiferro- 
magnets) rotates in the "easy-plane" and the propagation 
vector is perpendicular to this plane. In these noncen- 
trosymmetric magnets such spirals are stabilized by Lif- 
shitz invariants with gradients along the "hard-axis" . In 
our model, however, the Lifshitz invariants include only 
gradients in the basal plane. Correspondingly, the chiral 
modulations in these systems have propagation directions 
only in the basal planeOllj 

For easy-axis systems (K > 0) the critical surface 
D c (h±,h z ) has a more involved shape. In this case the 
lowest value of D c equals zero. This is reached at the 
spin-flop field (0, ±h ) (Figs. @(b), |(b)). Thus, near the 
spin-flop transition the modulated states arise at arbi- 
trarily small values of D. This unusual situation is due 
to the particular evolution of the potential profile ( |l6| ) 



(a) 



o 




(b) 
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FIG. 8: Contour lines D c (h±, ft, z )=const of the critical sur- 
faces for modulated states, fa) K < surface given by (pH[); 
(b) K > surface according to (||), @, and @. (The 
direction of increasing D c is indicated by dotted arrows.) 



in a magnetic field directed along the tetragonal axis 
(Fig. |^(b)). In this case the uniaxial anisotropy and the 
applied magnetic field have competing influence on the 
magnetic structure. While the easy-axis anisotropy ori- 
entates the staggered magnetization along the tetragonal 
axis, the applied field orientates it perpendicular to this 
axis. An increasing magnetic field in the region h < ho 
gradually decreases the potential barrier between the 
states of the antiferromagnetic phase with antiparallel 
staggered magnetization (Fig. ||(b)). When the spin-flop 
is reached the applied field completely cancels the influ- 
ence of uniaxial anisotropy, and the potential ( |l6|) equals 
zero for any orientation of the staggered magnetization. 
This infinite degeneracy of magnetic states is artificial 
because of the neglect of higher order anisotropy con- 
tributions. A fourth-order uniaxial anisotropy K-i sin 4 9 
removes this degeneracy. We may generally state that in 
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noncentrosymmetric easy-axis antiferromagnets near the 
spin-flop field the potential barrier between minima (and, 
therefore, the critical values of D) is determined by the 
much weaker fourth-order anisotropy constant Ki. We 
add two remarks here: (i) For centrosymmetric eas y-ax is 
antiferromagnets (models I.) and II.) of subsection [I A ), 
the situation near spin-flop transitions iSj.nhjz.sically com- 
parable and has been studied in detail.Ej'Efl There, the 
relation between interactions and homogeneous magnetic 
states is well understood, (ii) In cubic helimagnets,E3 
where magnetocrystalline anisotropy is represented only 
by fourth-order terms, no suppression of the modulated 
states has been observed, rather chiral modulations ex- 
ist in the complete region of existence of magnetically 
ordered states. Estimates based on the physical origin 
of these magnetic energy contributions yield so weak 
threshold values that modulated chiral states-Jiiear the 
spin-flop field should be expected generally!!!] There- 
fore, this competition between antisymmetric exchange 
and anisotropics makes the easy-axis noncentrosymmet- 
ric antiferromagnets particularly interesting systems for 
a search for and investigation of modulated chiral states. 



IV. LOCALIZED CHIRAL STRUCTURES 




DWII 



DWI 



FIG. 9: Structure of domain walls in the canted phase given 
by the turn angle A8 as function of distance from the wall- 
center £ = for different values of the parameter v. The 
180-degree domain walls of the AF phase (u — 0) are de- 
formed in the canted phases either into walls with a decreased 
value of A9(£ — > oo) (DW I, profiles with continuous lines) 
or an increased A#(£ — » oo) (DWII, profiles with dashed 
lines). DW II are transformed into 360-degree walls in the 
limit v = 1. 



In this section we consider the influence of the 
Dzyaloshinskii-Moriya interactions on localized magnetic 
defects within homogeneous magnetic configurations. As 
there is a wide variety of possible defect structures in 
the different phases, we will present only a few examples 
to demonstrate the general principles which rule defect 
structures for the noncentrosymmetric antiferromagnets 
under the influence of chiral couplings. A formal math- 
ematical description of isolated, planar one-dimensional 
defect structures, i.e. domain walls, was already devel- 
oped above in section IIIC. In the next subsection, we 



discuss physical importance and properties of such do- 
main walls for antiferromagnetic systems described by 
our model. Subsection IVB is devoted to linear two- 
dimensional defect structures, i.e. vortices. 



A. Domain walls or kinks 

Planar defects (domain walls or kinks) are commonly 
observed magnetic localized, states in many classes of an- 
tiferromagnetic materials.ES They separate homogeneous 
states with different degenerate directions of the stag- 
gered magnetization. An example is provided by 180- 
degree domain walls between regions with antiparallel 
staggered magnetization, i.e. different antiferromagnetic 
phases, in easy-axis antiferromagnets. In the antiferro- 
magnets under discussion, the rotation of the vector n 
within a domain wall is accompanied by oscillation of the 
total magnetization. The spin arrangement in such do- 
main walls is similar to that in the corresponding spirals 
(Fig. |^). Rotation of n as in a Bloch wall with longi- 



tudinal modulation of the vector m (Fig. H(a)) should 
occur in noncentrosymmetric antiferromagnets belong- 
ing to crystallographic class D^d- Neel- wall- like struc- 
tures with transversal oscillation of the magnetization 
correspond to the noncentrosymmetric antiferromagnets 
from class C\ v (Fig. ||(b)). The inhomogeneous chiral 
Dzyaloshinskii-Moriya interactions do not influence the 
structure of the domain walls (see above, remarks fol- 
lowing eq.(p0[)). However, the domain wall energies do 
depend on the rotation sense according to (|i|). As dis- 
cussed above, the modulated structures have a fixed sense 
of rotation that corresponds to a decrease of total energy 
compared to the homogeneous states. Spirals with op- 
posite sense of rotation are unstable (even with respect 
to the homogeneous states) and never arise in real sys- 
tems. Contrary to this, domain walls with disadvanta- 
geous sense of rotation, although increasing the energy, 
should be found within these antiferromagnets with sim- 
ilar probability because domain walls in antiferromag- 
nets have mostly "kinetic" origin in contrast to domain 
structures in ferromagnets, i.e. antiferromagnetic do- 
main structures are formed during the transition to the 
ordered states or as a result of reorientation transitions. 
These processes are largely independent of domain wall 
energies. 

The structure of such domain walls can be derived by 
integration in (p2[). Here, we restrict ourselves to one 
example of a practical calculation. We obtain the struc- 
ture and characteristic parameters of the domain walls 
in the canted phase (see Figs. |l|, |J). Within all regions 
of their existence {y < 1) the equilibrium states are sep- 
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FIG. 10: The domain wall widths A as functions of the 
parameter v for the two wall types in the canted phase, (a) 
DW I (b) DW II, according to (§§). 
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arated by two types of barriers in the potential profile 
$(#) (Fig. ||(a)). Correspondingly there exist two types 
of domain walls in the canted phase. The first low energy 
domain wall (DW I) separates homogeneous states with 
9 — arcsin z^ and 9 = it — arcsin z/; and DW II correspond- 
ing to the higher potential barrier separates states with 
9 = arcsin z; and 9 = —ir — arcsin z/ (Fig. ||(a)). Evaluat- 
ing the integral (^2|) with <I> from ([!(]) yields the following 
results 



sin 9 



z/cosh(£Vl - v 2 ) ± 1 



3h(£vT 



± v 



(28) 



Wall structures for both domain wall types in the canted 
phase with varying v are displayed in Fig. ^. The effec- 
tive thickness of domain walls A is usually determined 
as a distance between points where the tangent at the. 
inflection point intersects the lines 9 — 9\ and 9 = 92x3 
For our example this definition yields the following ex- 
pression 



A = '^"f^„ 



7r =p 2 arcsin v 



(29) 



(In fl2q ) and (j||) the upper/lower signs correspond to 
DW I/DW II.) The dependence of the wall thickness 
A on v for both types of walls is shown in Fig. [l(]. 
For increasing v the difference between magnetic con- 
figurations in the adjacent domains separated by DW I 
(A9 = 7r — 2 arcsin v) and the potential barrier (APhi = 
$(7r/2) — $(#i) ~ (1 — v) 2 ) gradually decreases while the 
thickness of the wall increases. At the critical point of 
the transition into the weak ferromagnetic phase, v = 1, 
the difference between magnetic states in the domains 
disappears and the wall spreads out without bounds. 
For DW II the potential barrier A<I> ~ (1 + v) 2 ) and 
A9 = tt + 2 arcsin v increases with increasing v. At the 
critical point, v = 1, these walls transform into 360- 
degree domain walls (Figs. || and |I(i| ). 
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FIG. 11: Vortex structure for antiferromagnets with C„„ 
symmetry (a), (c),(e) and D2d symmetry (b), (d), (f). (a) and 
(b) distributions of staggered vector 1. (c) - (f) projections 
of 1 and oscillating total magnetization m in the basal plane, 
respectively. 



Structures and parameters for 180-degree domain walls 
in the antiferromagnetic and spin-flop phases can be de- 
rived in a similar way. All these domain walls may play 
the role of nucleation centers during the transition from 
the homogeneous to modulated states. On the other 
hand, as demonstrated in the previous section, at a tran- 
sition into the homogeneous state the spiral states break 
down into a system of isolated plane walls. 



B. Vortices or skyrmions 

Linear magnetic defects are another type of topological 
excitations that can exist in noncentrosymmetric mag- 
netic crystals due to the stabilizing effect ofpthe inho- 
mogeneous Dzyaloshinskii-Moriya interactions 113 At zero 
magnetic field and in fields applied along the tetrago- 
nal axis the model ([|) is invariant to rotation about the 
z-axis. Solutions for the vector n(r) axially symmet- 
ric in the basal plane and uniform along the tetragonal 
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FIG. 12: Delocalized vortex structure in the spin-flop phase 
in crystals with C„ v symmetry. Staggered vector 1 is shown. 



axis, i.e. vortices, obey this symmetry. As an example 
of such localized states we consider an isolated vortex 
in the antiferromagnetic phase (K > 0, h is parallel to 
the tetragonal axis and smaller than the spin-flop field 
h < ho). We assume that the staggered magnetization 
is oriented parallel to the z-axis on the vortex axis and 
rotates into the antiparallel orientation with increasing 
radial distance from the vortex core. It is convenient 
to introduce cylindrical coordinates for the spatial vari- 
ables, r = :ro (p cos (/?, psinip, z), in the expression for the 
energy (|o|). (As earlier in the case of the spirals, we use 
length units x = y/AX/]K\.) The analysis of the energy 
functional ( |Io|) shows that the problem has solutions 9(p) 
with 

tp = <p for class C nv , i> — ir/2 — ip for class Did- (30) 



For C m symmetry the solution (|30j) describes a cycloid- 
like rotation of the staggered magnetization vector 1 
(Figs. 0(a) and (c)). In the case of Did symmetry the 
vortex has more a sophisticated structure (Figs. pT|(b) 
and (d)). The rotation of the staggered magnetization 
in the vortices is accompanied by in-plane oscillations of 
the total magnetization (Figs. [0(e), (f)) as described by 
eq.fl). 

The equilibrium distribution 0(p) is determined from 
the differential equation common for both classes 



dp^ 



1 d6 
pdp 



sm v cos f 



(31) 



AD sin 2 9 



nD, 



o P 



h 2 



sin 9 cos 9 = 



with boundary conditions 9(0) = and 9(oo) = tt for 
localized vortices in the antiferromagnetic phase (h < ho, 
Fig. [Tl| ), or 9(oo) = ir/2 for delocalized vortices in the 
spin- flop phase (h > ho, Fig. |l2|). These different bound- 
ary conditions result in an important physical difference 
between these two cases of antiferromagnetic vortices: 
For the localized vortices the homogeneous equilibrium 
state is established everywhere for p — > oo and the in- 
homogeneity is localized in the vortex core. In the case 




FIG. 13: Vortex profiles 9(p): solutions of (|3lj) for different 
values of D / Do in zero applied field. 



of the spin-flop phase, the vortex structure at p — > oo 
is inhomogeneous with 9 = tt/2, but the angle ip rotates 
through the a full circle from to 2ir. Therefore, these 
vortices are named delocalized. Eq. (|3l| ) has solutions 
only when D is smaller than the critical values for the 
transition to the modulated phase D c (h) given by j2rj). 
Typical solutions for 9(p) are plotted in Fig. As 
D approaches the critical value D c the vortex expands 
without bounds. Eq. (|3l]) functionally coincides with 
the equations for isolated vortices in other models with 
Lifshitz invariaftte, i.e. models for noncentrosymmetifie 
ferromagnetsJl3e3 for other classes p£antiferromagnets,E2l 
as well as for chiral liquid crystals.^ For detailed analy- 
sis of eq. ([Il]) and discussion of the related questions see 
these papers. 

The vortices or skyrmions considered here (Figs, [ll], 
|l~2| , |l~3| ) are non-singular linear defect structures. They 
belong to topological defects studied in many fields of the 
modern physics—Similar topological objects arise in su- 
perfluid helium,ciip two-dimensional electronic systems 
(Hall skyrmions) p3 or in nanomagnetic materials. c3 It is 
important to mention that there is a fundamental cor- 
respondence between these theoretical models.c^l In the 
isotropic case (wd = w = in ([To|)) the equation for 
the vortex has analytical solutions which are well-knows, 
as Belavin-Polyakov-solutions for nonlinear i7-models.c£l 
These solutions turned out to be unstable in centrosym- 
metric magnetic crystals and collapse spontaneously un- 
der the influence of anisotropic internal interactions or 
applied magnetic fields. Thus, the Lifshitz invariants 
are crucial for stabilizing these vortes structures in non- 
centrosymmetric magnetic crystals. c3 Hence, such low- 
symmetry magnetic crystals are interesting and impor- 
tant systems for investigations of general properties of 
vortices. 

In the spin-flop phase the vortex states have delocal- 
ized character (Fig. [12]). They are similar to vortex states, 
in liquid helium or some textures in liquid crystals .eZI 
Such vortices for noncentrosymmetric antiferromagnets 
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with d = have been investigated in Ref. [50| They read- 
ily form localized vortex pairs similar to those responsible 
for Berezinskii-Kosterlitz-Thouless transitions.^ An ap- 
plied magnetic field deviating from the tetragonal axis vi- 
olates the axial symmetry of the system. Then, the two- 
dimensional localized states are expected to have vari- 
ous elongated shapes similar to those observed in chiral 
liquid crystals .c3 The Lifshitz invariants can also stabi- 
lize three-dimensional localized states (as free spherulites 
or drops). Ej Up to now no experimental observations or 
theoretical investigations of such structures have been re- 
ported. 

Concluding this section we draw attention to an impor- 
tant difference between the localized states in our model 
and those in other magnetic systems. In noncentrosym- 
metric antiferromagnets with weak-ferromagnetism due 
to the oscillating weak magnetization in the basal plane 
the domain walls and the vortices are susceptible to the 
influence of applied magnetic fields. For this reason, non- 
centrosymmetric antiferromagnets are convenient for the 
study of phase transformations and the dynamics of such 
nonlinear localized excitations. 



V. RELATION TO EXPERIMENTAL 
OBSERVATIONS 

The known noncentrosymmetric antiferromagnets with 
weak ferromagnetism include two tetragonal crystals. 
First, for Ba2CuGe207 belonging to the crystallographic 
class 42m (-DmJ, chiral modulations were discovered 
five year ago.E3 For the crystallographic class 4mm 
(Ci v ) only the antiferromagnet K2V3O8 has apparently 
been investigated and no modulated states have been 
found yet l 18 li 5C l Here, we shortly review experimental data 
for these two antiferromagnetic compounds within the 
framework of our theory. In the last subsection we com- 
ment on related experiments on noncentrosymmetric an- 
tiferromagnetic crystals. 



A. Ba 2 CuGe 2 7 

A chiral spiral with propagation vector in the basal 
plane and period length of about 37 unit cells was found 
as magnetic ground state for Ba2CuGe207 (space group 
P42im).E3 It was also found that a rather strong mag- 
netic field applied along the tetragonal axis induces a 
transition into a homogeneous state.ED The field depen- 
dencies of the period and the magnetization reported 
are in quantitative agreement with theoretical results of 
Ref. |. It appears that there is no local minimum of 
the period length in dependence on the strength of a 
magnetic field applied along the tetragonal axis. This 
implies that the uniaxial anisotropy of this crystal is of 
easy-plane type (K < 0) (cf. Fig. ^). To analyze their 
experimental data the authors explored the model with 
d = 0. In their experiments there is no indication of 



effects related to weak ferromagnetism. 

Further detailed investigations in magnetic fields ap- 
plied along other directions are required to determine 
the character of the uniaxial anisotropy and the values 
of the other characteristic parameters of the magnetic 
system within the general phcnomcnological expression 
(||) for the energy. 

B. K 2 V 3 8 

For this compound, at a temperature of 2 K (the 
Neel temperature is about 4 K) the magnetization 
curves in a magnetic field along tetragonal axis and. 
in the basal plane indicate reorientation transitions .Ea 
These transitions are similar to those earlier observed in 
ccntrosymmetric antifewpjnagnets with weak ferromag- 
netism, e.g. hematite£-3c3 The authors conclude from 
neutron diffraction experim ent .that there are no indica- 
tions of chiral modulations.t§E2l According to the results 
of our theory such a situation may take place for easy-axis 
systems (K > 0) with weak chiral interactions (D < Do). 
As was discussed above, in the vicinity of the spin- 
flop field the criterion for the stabilization of the modu- 
lated states is considerably weakened. Thus, the search 
for modulated states in this system should be started 
from thorough investigation near the spin-flop field. We 
add that there are two other similar noncentrosymmet- 
ric vanadium oxides Rb2V30s and (NH^VsOs which 
are suaposed to possess antiferromagnetic order below 
10 K0 They could be investigated in search for effects 
of chiral interactions. 



C. Other noncentrosymmetric antiferromagnets 

The copper metaborate CuB 2 04 (space group I42d 
(£>2^)) belongs to the same noncentrosyppetcjc class 
as Ba2CuGe207. However, according totiaoO it has 
a more sophisticated four-sublattice antiferromagnetic 
structure with in-plane anisotropy. A long-periodic mod- 
ulated state has been observed in this crystal for a cer- 
tain temperature rangeSEil Finally we mention here two 
other noncentrosymmetric antiferromagnets: a modu- 
lated chiral state has been observed in the noncentrpsym- 
metric antiferromagnet BiFeC>3 (space group i?3c)E3For 
CuFeS2 (space group 742d (D^)) antiferromagnetic or- 
der was reported to exist, however, no details about the 
magnetic structures are given E3 

VI. CONCLUSIONS 

In this paper we shepstthat the noveLantiferromagnetic 
crystals Ba2CuGe 2 07ji3 and K2V3O8E3 in spite of the re- 
ported difference in their magnetic properties belong to 
a common, previously unknown class of magnetic crys- 
tals: noncentrosymmetric antiferromagnets with weak 
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ferromagnetism. The phenomenological expression for 
the magnetic energy of such systems including all in- 
teractions allowed by symmetry (^) can be reduced to 
the functional (||) which describes the orientation of the 
staggered magnetization and can be considered as gen- 
eral model for two-sublattice antifcrromagnets. It in- 
cludes, as specific cases, all main classes of antiferromag- 
netic crystals (collinear antifcrromagnets, antifcrromag- 
nets with weak ferromagnetism, noncentrosymmetric an- 
tifcrromagnets without and with weak ferromagnetism). 
Further, by using realistic assumptions about the rela- 
tive strengths of the phenomenological constants in jlO| ) 
the problem has been reduced to the case that the ro- 
tation of the staggered magnetization is restricted to a 
certain fixed plane. This simplification yields a represen- 
tative and realistic approximative model replacing the 
general model @. It is amenable to a complete analy- 
sis of the possible solutions for magnetic structures. The 
boundaries of their existence in parameter space could 
be calculated in all detail and a clear physical picture 
of the formation and evolution of these magnetic states 
is achieved. Due to the unique combination of those in- 
teractions inducing weak ferromagnetism and those sta- 
bilizing modulated chiral states, a rich variety of new 
modulated and localized structures was found to exist in 
this class of magnetic crystals. In these inhomogeneous 
states chiral rotation of the staggered magnetization is 
always accompanied by oscillations of a weak magnetiza- 
tion component in the basal plane (Fig. ^). The modu- 
lated states in these systems can be realized as structures 
with the propagation vector along certain in-plane direc- 
tions (spirals). We remark that another type of solutions, 
two-dimensional modulated phases, so-called vortex lat- 
tices may also exist. In noncentrosymmtric ferromagnets 
they are thermodynamically stable .-wider applied fields 
in certain region of the phase spaceJla In Ref. ^ vortex- 
lattices in antifcrromagnets lacking inversion symmetry 
with d = 0, have been studied theoretically. Nucleation 
of such vortex lattices during the transition from the spin- 
flop phase is discussed in Ref. [3?]. However, it is still 
unknown whether these vortex lattices can be thermody- 
namically stable in antiferromagnetic materials. 

We have described one-dimensional localized struc- 
tures (domain walls or kinks) separating domains of 
homogeneous states. These differ from similar objects 
found in many other classes of magnetic materials by 
oscillations of the local net magnetization and the de- 
pendence of their energy on the sense of rotation for the 
staggered magnetization within the wall. These peculiar- 
ities of their properties should be accessible to experimen- 



tal verification. Furthermore, two-dimensional localized 
structures with finite sizes (as axisymmetric vortices in 
the antiferromagnetic phase) are possible topological de- 
fects in these systems. They are stabilized only due to 
the chiral interactions. 

In this paper we have deliberately avoided a detailed 
investigations of the full model (||). Instead, by intro- 
ducing a simplified model, we have described the general 
features of the magnetic properties in noncentrosymmet- 
ric tetragonal antifcrromagnets. We expect that this phe- 
nomenological description will provide a guide for further 
detailed experimental investigation of the known non- 
centrosymmetric tetragonal antifcrromagnets and for a 
search of new crystals belonging to this group. 

We also briefly indicate here possible further direc- 
tions of the theoretical investigations. Fourth-order 
anisotropics are needed to describe orientational pro- 
cesses in the basal plane and the peculiarities of mag- 
netic properties near the spin-flop field. Future theo- 
retical investigation also should include the stray field 
effects responsible for multidomain states near the first- 
order phase transitions. Similar investigations within the 
general model (||) pose a much more complex and chal- 
lenging task. This functional can be considered as a gen- 
eralized version of the nonlinear tr-model, one of the ba- 
sic models in the theory of nonlinear physics and soli- 
tons. It is rela ted jt O |. m aay other models in condensed 
matter physics. BEaaOea The further development of 
the theory should involve the investigation of vortices 
and vortex. 1 attices as done for other noncentrosymmet- 
ric models.li3E] Similar multidimensional localized solu- 
tions of nonlinear field equations are i nten sely studied in 
many other fields of modern physics rffi'TEj Nucleation 
and evolution of such one-dimensional and two dimen- 
sional modulated patterns have deep physical relations to. 
similar patterns in superconductivity,E3 liquid crystals ,EJ 
and other condensed .matter systems and even in modern 
cosmological models.ES 
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